Cyclotron resonance in antidot arrays 



O 
O 
(N 

< 



D. G. Polyakov,!'* F. Evers,i and I. V. Gornyii-^'* 
^ Institut fiir Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 
^Institut fiir Theorie der Kondensierten Materia, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 

We study the dynamical properties of an electron gas scattered by impenetrable antidots in the 
presence of a strong magnetic field. We find that the lineshape of the cyclotron resonance is very 
different from the Lorentzian and is not characterized by the Drude scattering rate. We show that 
the dissipative dynamical response of skipping orbits, Sc{uj), is broadened on a scale of the cyclotron 
frequency lJc and has a sharp dip oc \uj — ujc\. For small antidots, S'c(w) is strongly modulated with 
a period equal to uic and has sharp square-root singularities for a series of resonant frequencies. For 
large antidots, Sc{i^) has a hard gap a.t uj < Uc between two sharp peaks, associated respectively 
with edge states and free cyclotron orbits. 
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I. INTRODUCTION 



Progress in controlled fabrication of semiconductor 
nanostructures has revived interest in quasiclassical 
features of transport in a two-dimensional electron gas 
(2DEG). Since the Fermi wavelength of electrons in high- 
mobility heterostructures is usually small compared to 
the characteristic spatial scale of inhomogeneities, the 
transport properties of the 2DEG retain signatures of 
the underlying quasiclassical dynamics of the particles. 
In particular, transport in ballistic mesoscopic systems, 
where electrons are scattered specularly on the boundary 
of the system, has been investigated in terms of quasi- 
classical dynamics in considerable detail In antidot 
(AD) arrays, potential barriers around the ADs can also 
be viewed as hard discs of size ~ 10 — 10^ nm that re- 
flect electrons specularly. If AD arrays are periodic (for 
a review see [§4)), the quasiclassical character of elec- 
tron dynamics manifests itself in pronounced geometri- 
cal resonances, which are associated with the periodicity 
and lead, in particular, to commensurability peaks in the 
magnetoresistance. On the other hand, random AD ar- 
rays (for experimental work on dc transport see, e.g., 
) represent a remarkable disordered system in which 
the statistics of fluctuations of the random potential is 
strongly non-Gaussian (in contrast to the familiar case of 
smooth disorder in high-mobility heterostructures, where 
the random potential at a given point is a sum of contri- 
butions from many impurities). 

On the theoretical side, the recent interest in the quasi- 
classical dynamics of a 2DEG is to a large extent inspired 
by a variety of "non-Boltzmann" quasiclassical trans- 
port phenomena that occur in disordered systems with 
large-scale inhomogeneities. The term "non-Boltzmann" 
means that these phenomena, while being essentially 
classical, cannot be described by the Boltzmann kinetic 
equation (i.e., the collision- integral approximation is in- 
sufficient). They are due to correlations of scattering acts 
at the points where quasiclassical paths self-intersect, 
which gives rise to memory effects, not captured by 
Boltzmann transport theory. Most noticeably, the non- 
Markovian kinetics yields a wealth of anomalous mag- 



netotransport phenomena in low magnetic fields [p[-p|, 
induces adiabatic localization of electrons in strong fields 
[|l2|,Q, and leads to a peculiar behavior of the magne- 
torcsistivity in the Lorentz model of hard-disc scatter- 
ers It is important that the quasiclassical non- 

Boltzmann corrections dominate over the quantum ones 
in systems with long-range disorder. 

The non-Markovian character of kinetics that leads to 
the anomalous dc transport manifests itself also in the 
cyclotron resonance (CR). In particular, the adiabatic 
localization ||l2|,Q is predicted to yield a peculiar 
shape of the CR line: a narrow peak related to chaotic 
dynamics of delocalized electrons on top of an inhomo- 
geneously broadened background coming from adiabati- 
cally localized electrons. The latter contribution domi- 
nates for large B and gives the linewidth which depends 
nonmonotonically on B. This intricate picture should be 
contrasted with the case of white-noise disorder, where 
the linewidth is given by the scattering rate ||l^ . 

In this paper, we consider the CR in AD arrays. We 
demonstrate that the CR lineshape is very much different 
from the Lorentzian suggested by the Drude theory. The 
peculiarity of the dynamical response of the AD system 
is related to two factors which become important with 
increasing B: formation of "skipping orbits" of electrons 
bound to ADs and suppression of scattering for other 
electrons that do not participate in the process of skip- 
ping. We show that the CR is not characterized simply by 
the Drude scattering rate. Specifically, cyclotron orbits 
not colliding with ADs yield an infinitely sharp CR line, 
whereas skipping orbits give a contribution broadened 
on a scale of the cyclotron frequency cUc- The skipping- 
orbit contribution exhibits a remarkably rich behavior as 
a function of frequency w: it has a sharp dip cx jtj — c^d 
at cj = Wc and, in the case of small ADs, is strongly mod- 
ulated with a period equal to uJc- The modulation yields 
exact zeros of the CR response in a dilute AD array. In 
addition to the zeros, a series of sharp singularities is de- 
veloped in the wings of the CR line. For large ADs, the 
dynamical response has a hard gap at w < Wc between 
two sharp peaks, associated respectively with edge states 
and free cyclotron orbits. 
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The paper is organized as follows. In Sec. ||, we derive 
an exact expression for the CR lineshape in the limit of 
large B, which is then analyzed in two essentially dif- 
ferent cases of small (Sec. and large (Sec. |V|) ADs. 
In Sec. we consider moderately strong magnetic fields 
(moderately in the sense that skipping orbits bound to 
different ADs can overlap). We add remarks bearing on 
experiment and the role of electron-electron interactions 
in Sec. VI. Throughout the paper, the analytical calcula- 
tion is complemented by results of numerical simulations. 



II. LARGE-B LIMIT: INSULATING PHASE 

We start by considering the limit of large B, namely 
ngR^ <C 1, where ns is the concentration of ADs, Rc 
the cyclotron radius. We model ADs by hard discs and 
assume that the array of ADs is dilute, i.e., nga^ ^ 1, 
where a is the radius of the discs. In this model (known 
as the Lorentz gas), if Rc is smaller than a critical value 



of the order of n 



-1/2 



all quasiclassical trajectories get 



localized and the dissipative dc conductivity axx = 0) 
vanishes at zero temperature exactly |l^Jl^. The lo- 
calization is developed through the formation of discon- 
nected clusters of trajectories that do not extend beyond 
a finite area. In random AD arrays, the metal-insulator 
transition is "second order" , so that as one goes deeper 
in the insulating phase with increasing B, the size of the 
critical clusters decreases continuously. Eventually, far 
away from the critical point, for nsR^ <^ 1, trajectories 
that collide with two or more ADs become very rare. In 
this limit, most electrons do not collide with ADs what- 
soever and give a (5-function CR line at uj = ojc- Most of 
those that collide move in skipping orbits around a single 
AD. It follows that for large B the dynamical response 
at (jj ^ uJc is determined by the skipping orbits. Clearly, 
this conclusion is true both for random and periodic AD 
arrays. 

Let us calculate the dynamical response associated 
with skipping orbits. Since in phase space of the Lorentz 
gas there is a well-defined separatrix between free cy- 
clotron orbits and trajectories colliding with ADs, we 
write the dissipative conductivity as a sum of two terms, 

Rea^A^) = ^e^Po[pDf{cu) + {1 ~ p)Dc{uj)] , (1) 



where 



p = exp(-27r/tJc'ro) 



(2) 



is the probability to close the cyclotron orbit without suf- 
fering a collision, tq = l/2vpnsa the collision time, po 
the density of states for free electrons at _B = 0. The sym- 
metrized functions Df^d'jj) = Vp[Sf^ci^) + '5'/,c(~w)]/2 
are the velocity-velocity correlators for free electrons 
and electrons colliding with ADs, respectively; vp is the 
Fermi velocity. In the above, we have assumed that many 



Landau levels are occupied, so that kpRc ^ 1, where kp 
is the Fermi wavevector, and that ADs are large enough, 
in the sense that kpa ^ 1, which means that ADs scat- 
ter electrons specularly. It is also important to us that 
under these conditions the density of states of electrons 
scattered by ADs is not affected by the Landau quanti- 
zation and is given simply by (1 —p)po. 

The dynamical response to a circularly polarized (CR- 
active) perturbation is given by S'/,c(w). For free elec- 
trons we have >S'/(a;) = ttS{uj — lUc), while for electrons 
skipping around an AD 



Sci^jj) = I dt (cos[([x'c — 



20^ e(i-|-Ai-mT)])At,r , 



(3) 



where Q{t) is the step function and 







At 



27r 



T(r) 

dAt 



Or 



2Rca J\R^^a 



rdr (4) 



denotes averaging over the initial phase Ai and the dis- 
tance r between the centers of the AD and the cyclotron 
orbit. The normalization of the integral takes care 
of the "exclusion volume" free of electrons because of the 
presence of ADs. The angle of incidence 9{r) (the angle 
between the trajectory and the tangent to the surface of 
the AD at the collision point) and the time T(r) between 
two successive collisions read 



COS! 



cos ■ 



(r) = 
uJcT{r) 



r'-Rt- a' 

2Rca 
a'^ -Rl- 

2R~^ 



(5) 
(6) 



Since r is an integral of motion for orbits skipping around 
a disc, 6{r) and T{r) are the same for each collision. 

Doing integrals over t and At and summing over m in 
Eq. (^, we get an identically zero response at the cy- 
clotron frequency for colliding electrons, Sc{(jJc) = 0, and 



n— — oo 
,2 , 



r„sin 6'(r„) 



E l n »n 



2Rca{LOc-ojf ^ |/:(r„)| 



(8) 



for LO ^ LOc- Here /[^(r) = ^{uJc'-^)T{r) + 9{r) and r„((jj) 
are roots of the equation fi^{r)+Tm = 0. The summation 
in Eq. (^) runs over r„ that satisfy |i?c — a| < < Rc + a. 

Equation (0) tells us that the dynamical response is 
due to resonant orbits. There is no broadening of the 
contribution of each of the orbits and, accordingly, the 
shape of the CR line is given by the density of states 
of the resonant orbits. The meaning of the resonance 
condition u) = ijjc + 2(0 -I- ■nn)T~^ is that the change of 
the total phase (cyclotron phase + phase of the ac field 
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+ scattering phase) should be a muhiple of 27r between 
two coUisions. It is worth noting that, contrary to what 
one might expect, the resonant orbits are not periodic 
two-dimensional orbits. While periodic orbits and tra- 
jectories close to them are indeed important in periodic 
arrays at weak B, where they give rise 20| to geomet- 
ric resonances, they do not play any particular role in 
Eqs. (0),(||)- This conclusion should be contrasted with 
the approach of Ref . |^ , where the "periodic-orbit anal- 
ysis" iQ, apphcable to the case of weak B, was extended 
to describe also the experimentally observed dependence 
of the position of the CR line on B in an AD array for 
skipping orbits bound to a single AD. Our calculation 
shows that this procedure is not justified. 

To analyze Sc{uj), it is convenient to rewrite Eq. (|^) 
and the resonance condition in terms of the function 
9lj{S) = fuj[r(0)] defined on the interval < 9 < tt: 



2{lUc - Lu)'- 



EhUi On 



g^i9) = -{Lu,~Lo)Ti9)+9 ; 



(9) 

(10) 
(11) 



The time between collisions T{9) behaves in an essen- 
tially different way depending on whether the ratio Rc/a 
is larger or smaller than unity. 



III. SMALL ANTIDOTS 



Let us first study the case i?^ ^ a by expanding T{9) 
m powers oia/Rc. Combining Eqs. (|),(|) we get 



T{9) 



27r 



2a 



sin 9 sin 29 

2Rr 



(12) 



with T(0) = T(7r) = 2tt/uJc- To find the roots 6'„ at 
|w — cjcl vp/a, one can retain only the first (unper- 
turbed) term in Eq. ([T^), which yields, for a given uj, a 
single solution 9n — T^iyij^c — [^^/^c]) and glj{9n) — 1- 
Here [lu/lOc] is the integer part of uj/lUc- Substituting 
these expressions in Eq. (B) we find 



S'c(w) = 



2{lJc-uj)^ 

2 I , ,2 



sm 



^ t:{lu — nujc) 



Dc{u)) = Vp- ■ 



(13) 
(14) 



for nojc < uj < {n + l)u!c and \n\ <C Rc/a- We thus see 
that the CR line exhibits a strong modulation, namely 
Sc{oj) has zeros at w = noJc- Note that the zeros are 
exact even if higher-order terms in powers of a/ Rc are 
taken into account. These zeros correspond to resonant 
orbits with 9 ^ Q and 9 ^ tt, which go along the tangent 
at the point where they touch the AD. In the vicinity of 
the zeros Sc(oj) behaves as \oj — nuJc\^ for any n ^ 1, in- 
cluding the dc limit n — 0. At the point uj — lOc-, Sc{lo) 



vanishes as jw — wd- The envelope of the oscillations ( |T^ ) 
falls off with increasing lu as (uj — ujc)~'^, similarly to the 
conventional Lorentzian. 

The total response of the 2DEG scattered by ADs of 
small radius a {Rc/a ^ 1 and |a;| ^ vp/a) is thus a 
sum of two parts: a sharp peak associated with Sf{uj) 
and a series of broader peaks given by Sdco), whose 
width is ~ LOc- The oscillatory behavior of Dc{lo) is 
shown in Fig. |^. One can see that the first two peaks 
of Dc{uj) that occur between w = and lj = 2u)c are 
much higher than those for larger w. As a result, the 
response of skipping orbits with Rc/a ^ 1 looks like a 
double peak split up at = Wc- Note that the resonant 
value of 9 tends to tt if one approaches uJc from the left 
{lj LOc — Q), whereas it tends to if one does so from 
the right {lo — > + 0). Accordingly, the two parts of 
the double peak for 1 < lj/lUc < 3/2 (where the reso- 
nant angles < 9 < tt/2) and 1/2 < lo/lUc < 1 (where 
7r/2 < < tt) differ in the direction in which the guiding 
centers of the resonant skipping orbits rotate around the 
AD (see the insets in Fig. 1). 
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FIG. 1. Dynamical response of skipping orbits for 
Rc/a ^ 1 and uj <^ vp/a exhibits oscillatory behavior with 
a characteristic double peak around zero ai uj — cuc- Dashed 
line: Dc{uj) in units of vp/ujc according to Eqs. (p^,(p^; 
solid line: the numerical simulation for Rc/a ~ 11.3. The 
insets illustrate the different sense of rotation for resonant 
skipping orbits with 1 < uj/uJc < 3/2 (upper panel) and 
1/2 < uj/ujc < i (lower panel). 

The wing of the principal double peak [Fig. |^ exhibits 
nontrivial behavior. As lu increases, Sc{uj) in intervals be- 
tween two adjacent zeros gets more and more asymmet- 
ric with maxima shifting towards the higher-|a;| bound- 
ary of the intervals. Eventually, when \uj\ reaches the 
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critical frequency vp/a, the behavior of S'c(u;) acquires 
quaUtatively new features. At jw| > vp/a, gui(0) as a 
function of 9 becomes nonmonotonic, which can be seen 
from Eqs. (|lo|) , ([l^) . This leads to the appearance of mul- 
tiple roots 9n for a given oj (in contrast to a single root 
at |a;| < vp/a)- Because of the multiple roots, the zeros 
in Sc{w) disappear for |ti;| > vp/a, since the condition 
= or TT now cannot be met for all roots simultane- 
ously (so that there is a finite number of zeros, namely 
2Rc/a zeros at Rc ^ a). Moreover, the nonmonotonic 
behavior of g^^ (9) yields singularities in the lineshape as- 
sociated with resonant orbits for which the derivative 
g'^{9n) vanishes in the denominator of Eq. (0). 




FIG. 2. Dynamical response of skipping orbits with 
Rc/a ^ 1 in the tail of the principal double peak (shown 
in Fig. The oscillating curve represents the product 

Dcioj) X (lli/ujc)^ in units of vp/uJc as obtained from the nu- 
merical simulation for Rc/a ~ 11.3. The inset magnifies the 
region near uj/ujc = Rc/a, where zeros of Dc{uj) disappear 
and simultaneously a series of square-root singularities starts 
[Eqs. (^5|),(|l^)]. It also shows the resonant orbit correspond- 
ing to the singular frequency. 

Let us analyze the high-w limit \uj\ 3> vp/a. In this 
case, one can neglect the last term in the right-hand 
side of Eq. (|l^) and represent the equation for 9n at 
n ^ 1 in the form (uja/vp) siii9n + 7r(n — lo/uJc) = 
with g'^{9n) — (uja/vp) cos9n. One sees that now for a 
given to there are 2\Lo\a/TTVp ^ 1 roots 9n- Transforming 
to the continuous limit in the summation over 0„ [Eq. (|^)] 
we obtain a regular part of Sc{lo) for 3> vp/a: 



-^c (^) ^ o 2 



(15) 



The most prominent feature on top of the smoothly vary- 
ing background ( p^ is the appearance at \uj\ > vp/a of 
sharp singularities ("spikes") in Sc{(^)- The spikes oc- 
cur every time g[j{9) vanishes for one of the roots 0„. 



Note that the singularity frequencies coincide for Sc{(^) 
and S'c(— w). At \uj\ ^ vp/a, gu){9) reaches maximum at 
9 = 9max — 7r/2 + a/Rc + vp/cua, close to 7r/2. It follows 
that the skipping orbits that yield the spikes hit the sur- 
face of an AD at almost a right angle, so that the center 
of the orbits flicks from one side of the AD to the other 
after each collision. Expanding around the maximum 
and substituting two (almost degenerate) roots 9n{Lu) of 
Eq. ( pi] ) which are close to 9 max, we find Sc{uj) in the 
vicinity of the singularities for ^ vp/a: 



Qsinq ( , ,\ / 



Vp 



up' \27ra|a;| 

1/2 



1/2 



6[(cj„ - Lo) sgno;] 



(16) 



The frequencies a;„ for |w| 3> wf/o and Rc^ a are given 
by ujn — Wc(n -I- 3/2)(l -I- a/nRc). Note that the period 
of the sequence of spikes is larger than (but close to 
ujc at Rc S> a). The behavior of S'c(tf) for large ua/vp 
is illustrated in Fig. ^. The square-root singularities at 
Lu — > ujn that appear for uj > vp /a (i.e., for uj /u)c > Rc/o-) 
are clearly seen. The region of a; ~ vp/a is blown up in 
the inset to show the sudden start of the series of singu- 
larities when Lo becomes larger than vp/a. 



IV. LARGE ANTIDOTS 

We now turn to the case Rc ^ a. As we will see, Sc{oj) 
shows completely different behavior in the two limits of 
large and small Rc/a. The function T{9) expanded in 
powers oi Rc/a reads: 



T{9) 



2{tt - 61) 2Rc 



ujca 



iRc 

4a 



sin 29 



(17) 



with r(0) = 2ti/lOc and T{ti) = 0. Clearly, if one sends 
Rc/a —9- 0, Sc{>jj) describes the dynamical response of 
trajectories skipping along a straight hue ("edge states"). 
Using ( p7| ) in this hmit, we get 1 -I- [Iwl/cjJ roots 0„ = 
■K — {n -\- 1)t:uJc/\lo\ with < n -f- 1 <[|a;|/a;c], for which 
9'u{(^n) = '^/'^c- Substitution in Eq. (O) yields 



2|a;|(wc 



[|w|/(^e]- 

E 

n=0 



sin^ ^ , ; — - (18) 



and Dc{uP) which is related to S'c(w) by Eq. (|l4|). Since 
the limits w — s- and Rc/a — > do not commute with 
each other, one should be careful about the behavior of 
Sc{uS) at ^ 0. In fact, 5'c(w) for the edge states has 
a (5(Li;)-term in the dc limit, so that Eq. ( p^ correctly 
describes the response of a large AD with Rc/a only 
at a; 7^ (see below). 

As follows from Eq. (|l8|), in a striking difference with 
the case of large Rc/a, there appears a gap with Sc(uj) — 
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for < |a;| < lOc- No dissipation occurs with increasing 
oj until exceeds loc- At the edge of the gap, Sc{io) 
vanishes hnearly as jwl — Also, in contrast to the case 
of large RJa, Sc{^) given by Eq. has neither zeros 
[cf. Eq. (p^J nor singularities [cf. Eq. (p^j for > uic- 
In the limit of large \u!\ ^ uic, replacing the summation 



in Eq. ( p^ ) by integration leads to Sc{u 
behaves according to Eq. ([l5|). 



which 




The total weight of the peak at small to is given by 
Jo^^'diJ Sc{uj) = J^d0 (20)-i sin3 ~ 0.485. 

The dynamical response of the 2DEG in the case of rel- 
atively large ADs with Rc/a ^ 1 is thus a sum of three 
peaks: two sharp peaks, one centered at low frequency 
uj ~ vp/a, the other at higher frequency uj = lUc (the 
latter is related to free cyclotron orbits), plus a broad 
peak of width ~ lOc whose edge coincides with lOc- The 
behavior of the contribution of skipping orbits, S'c(w), is 
illustrated in Fig. |[ 

Note that the hard gap in S'c(w) gets narrower but 
remains exact at small Rc/a. In fact, the hard gap sur- 
vives with increasing i?/a up until Rc becomes equal to 
a, at which point the geometry of skipping changes in a 
qualitative way. Specifically, &i Rc < a skipping orbits 
propagate in only one direction around the AD, whereas 
for Rc > a they propagate in both. As a result, the hard 
gap at < ujc disappears (transforms into a soft gap 
oc jw— Wcl) for Rc > a and simultaneously infinitely many 
singularities pop up at |a;| > uJc- With further increas- 
ing Rc/a, zeros in Sc{uj) proliferate in the finite range of 
frequency < vp/a, as explained above. 



V. MODERATELY STRONG B: 
PHASE 



METALLIC 



FIG. 3. Dynamical response of skipping orbits with 
Rc/ct <C 1 (edge states). The sharp low-frequency peak is 
separated from the broad peak at tj ~ tJc by a hard gap. The 
curve for uj > cuc is Sc{io) in units of uj~^ calculated accord- 
ing to Eq. dlq). The peak at < lo < t;j?/a is a solution 
of Eqs. (jl|)!@) for Rc/a = 0.1 (solid line) and Rc/a = 0.2 
(dotted line). 

As mentioned above, in addition to Eq. ( [T8| ) there is a 
S{lo) peak in Sc{uj) for edge states with Rc/a 0. To 
calculate the dynamical response of ADs of a large but fi- 
nite radius, one should take into account terms of higher 
order in Rc/a ^ 1 in the expansion (pT|), which leads to 
a broadening of the peak in a finite range of frequency, 
namely < w < vp/a. This broadening is governed by 
the root which obeys the equation 



oja 

VF 



sm {/-i 



(19) 



1 /g) , where we retained 



with (7^(6'_i) — uj'^^{uj — vf COS. 

only the term linear in Rc/a in Eq. (|l7|). We see that a 
solution of ( p^ exists only in the abovementioned inter- 
val of LJ. Within this interval 5'c(ti;) reads 



Sc{l~o) = 7, 



2vp cosO^i + Lua/vp 



(20) 



Combining Eqs. (|19[),(|20|) we obtain the asymptotic be- 
havior of Sc{lli) near the ends of the interval: Sdco) van- 
ishes as at ^ and as {vp/a — uiY^^ at ^ vp/a. 



The above analysis of the dynamical response of a sin- 
gle AD applies directly to the case of a strong mag- 
netic field, nsRc <IC 1- In this limit, a skipping or- 
bit is bound to a single AD for an infinitely long time. 
With decreasing B, the scattering processes that in- 
volve collisions of a skipping orbit with many ADs be- 
come essential and at a critical value of the parameter 
nsR^ 1 a metal-insulator transition occurs. On the 
conducting side of the transition, the skipping orbits be- 
come delocalized by hopping from one AD to another. 
Deep in the metallic phase, for ngi?^ ^ 1, the char- 
acteristic hopping rate is Tq^ and the dc conductivity 
axx{0) is given by the second term in Eq. ([!]) with |l5| ] 



DciO) 



-{x)/[1 + uj'^t'^{x)] parametrized by x ~ lOcTo, 



where t(oo) = tq and t(0) — 3to/4. At zero B, Boltz- 
mann theory works perfectly well in the hydrodynamic 
limit (ns — > oo, nga = const); but, due to the factor 
of 1 — p in Eq. (|l|), a^xiO) falls off as B^^ in the con- 
ducting phase, one power of B faster than in Boltzmann 
theory. One might think that, apart from this factor, 
the dynamical response at finite B will also be similar to 
that in the Boltzmann approach. The latter is simply the 
zero-B Lorentzian with a shifted frequency lu uj — tOc- 
In fact, however, the behavior of Sdto) in the conduct- 
ing phase is completely different from the Lorentzian, see 
below. 

While Sc{oj) in the dc limit shows the metal-insulator 
transition, at larger uj the function S'c(w) in the conduct- 
ing phase retains, provided uJcTq ^ 1, the main features 
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of the single- AD response with Rc/a ^ 1 [Eqs. (|T2|)- 
(|l6|)]. The oscillatory behavior of the lineshape with 
sharp dips at uj = nujc remains almost unchanged at 
'-^cTq 3> 1, since in this limit skipping orbits experience 
many collisions with a single AD before they "change" 
to another one. Clearly, in contrast to Eq. (|l3|), Sc{i^) 
has no exact zeros any more; however, the behavior of 
Sc{uj) is modified only in a close vicinity of the points 
Lu = ncjc- In particular, because of the broadening of the 
resonances due to the hopping between ADs, the linear 
vanishing of S'c(w) oc jcj — t^d near lu — lUc is cut off at 
|a; — o^d Tq^- On the other hand, the cubic vanish- 



ing of Sc{uj) 



matches Sc{0) ^ (w^tq) ^ in the 



dc limit, which establishes the scale ~ Tq^{lUcTo)'^^'^ 
on which the frequency dispersion of the conductivity 
becomes strong. A linear zero-frequency anomaly 122] 
Sc{uj) — Sc{0) cx \uj\ appears in the metallic phase. Notice 
that, because of the anomalously strong broadening of 
the CR line, the crossover to the regime (Tj.j.(a;) 3> crxxi^) 
occurs with increasing uj at much smaller (for uJcTq ^1) 
frequency than in the Drude regime, where the charac- 
teristic scale is cUc- 

The hopping between ADs also cuts off the square- 
root singularities (|l6|) [they only survive for isolated 
ADs, whose contribution is suppressed by the factor of 
exp(— 47rn5i?^)]. Note that the substitution ut uj+i/rQ 
in Eq. (16) is only correct ioi vp/a ut ^ vpRc/a^. 
For larger lo, the effective collision rate for the resonant 
orbits is renormalized, Tq^ — > Tq^, since after each two 
collisions with a given AD the center of the resonant orbit 
is shifted by a distance SR ~ RcVp/auj, which is much 
smaller than a for the large frequencies. It follows that 
To increases with lu: tq/to ~ a/SR, so that the singular- 
ities are cut off on the smaller scale of Tq^. Note that 
after the collision with another AD the center of the or- 
bit is shifted by ~ Rc{SR/ay^^ , which makes the orbit 
non-resonant. 

The smearing of the soft gaps around the zeros of Sd^j) 
can be clearly seen in the hydrodynamic limit 715 —> oo, 
a — > with nsa held fixed. In this limit, the kinetic prob- 
lem of finding axxi^^) allows for an exact solution p^ , p^ , 
which we reproduce for convenience in the following form: 



Sciuj) = 



1 



LUc 

X Re 



- U! 



1 — pe 



iT,{Lu) 



1-p 



U! — LUc 



with 



1 



(21) 



(22) 



Here g^{e) is given by Eq. (0) with T{e) = 27r/wc, 
since in the hydrodynamic limit Rc/a 00. Sending 
i^cTq — > 00 (i.e., p — > 1) generates a pole on the real axis 
of 6 in the integrand of the self-energy (22), which yields 
E(cj) ^ {TT/2To)sm^ Oe'" with 9 = 7r(w/wc - [^/wj)- The 



latter are precisely the 9's defining the resonant orbits 
in Eq. (|l^). Using the above expression for S(w) at 
^cTq ^ cxd in Eq. (21) indeed yields our Eq. (Esh, de- 
rived for the insulating phase. Clearly, no singularities 
[Eq. (p^j occur in the hydrodynamic limit. Note that the 
metal-insulator transition takes place in a dilute AD ar- 
ray at Rc/a ^ 1. Therefore, provided lUcTq ^ 1 (which 
is the conventional condition for a developed CR reso- 
nance), Sc{uj) is in fact described by Eq. ( p^ ) not only in 
the metallic phase but also in the critical region of the 
metal-insulator transition, except for a very close vicinity 
of the zeros of Sc{iu). 



Is/R,=28.28 R,/a=21.21 
14.14 42.42 

4.24 141.42 

2.12 282.84 




(o/co„ 



FIG. 4. Dynamical conductivity axxiLu) as a function of 
uj/ujc in units of the zero-B zero-tj value ao in a dilute antidot 
array as obtained by the numerical simulation. The transport 
mean free path at zero B is Is — 600a. Different curves cor- 
respond to different B. The ratio of Is/Rc changes from 2.12 
for the upper curve to 4.24, to 14.14, up to 28.28 for the lower 
curve. Even for the relatively weak field with Is/Rc — 2.12 a 
very sharp peak at lu = uuc is clearly seen. With increasing B, 
the broad double peak around lu = luc gets more pronounced, 
as does the oscillatory behavior of a^xiLu) for higher lu. For 
large Is/Rc = 28.28, the numerically obtained axxi^u) agrees 
well with Eqs. (§),(|l|),(0). 

In Fig. ^, we show results of our numerical simulation 
of ac transport in a very dilute AD array with nsa^ ~ 
0.62 • 10^"^, which corresponds to the mean free path 
Is = 3/8nsa ~ 600a, tuned through the metal-insulator 
transition by changing magnetic field. The upper curve 
corresponds to the smallest B, with Is/Rc — 2.12. One 
sees that, contrary to what one would expect from the 
conventional Boltzmann theory, although the dc value 
(^xxi^) is only ~ 6 times smaller than that at zero B, 
a very narrow peak at w = Wc is developed, associated 
with free cyclotron orbits. As B increases, this peak ac- 
quires more weight, while remaining very sharp. Next, in 
striking contrast to Boltzmann theory, the oscillations of 
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o'xxi'^), related to skipping orbits, become more and more 
pronounced with growing B. The double-peak structure 
around the linear gap near w = is clearly seen, as are 
deep minima at to equal to multiples of Uc- For large 
Is/Rc, the numerically obtained behavior of axxi^^) is in 
a good agreement with Eqs. (|l|) , , (|lj) . 

As follows from the comparison of Eqs. (13),(pi|),(22), 



the dissipation at cj ^ for finite uJcTq is due to resonant 
orbits ( |ll|) broadened by the hopping between different 
ADs. It is interesting to note that there appear two kinds 
of resonant behavior of the dynamical response in the 
metallic phase: first, the last factor in Eq. ( |2l| ) resembles 
a resonance at u; ~ broadened by ImE; on the other 
hand, the broadening itself is due to resonant orbits, for 
which the denominator in Eq. (p^) is close to zero. 



VI. DISCUSSION 

In the above, we have calculated the quasiclassical dy- 
namical conductivity of an AD array, completely neglect- 
ing electron-electron interactions inside the 2DEG. At 
zero B, in the case of Coulomb interaction in a normal 
metal, this is a well-controlled approximation governed 
by the large parameters kpag and kplg, where is the 
Bohr radius (and a characteristic screening length at zero 
B). At nonzero B, the situation appears to be much less 
trivial. Although in this paper we do not provide any 
treatment of the combined effect of the interactions and 
disorder on the CR line, a few comments are in order. 

Firstly, the Landau quantization enhances the role of 
the interaction. In particular, for weak B, the ground 
state of a disorder-free 2DEG is known [g3 24| to spon- 
taneously break translational symmetry in a partially 
filled (i.e. highest occupied) Landau level. It is elec- 
trons from this level and those from the highest fully 
occupied level that are excited in the CR at zero tem- 
perature. On the other hand, the CR at zero wavevector 
in a one-component system of particles with a parabolic 
dispersion without external inhomogeneities is insensi- 
tive to the electron-electron interaction (Kohn's theorem 
p5[). Hence, it is only because of the combined effect 
of disorder and interactions that the latter can affect the 
CR lineshape (e.g., |2^,^ and references therein). At 
present, it is unclear what is the resulting lineshape of 
the dynamical response for the system studied in [ p3|j24| 
(see also p8|-^) in the presence of disorder. However, 
since electrons skipping around hard-wall ADs do not 
experience the Landau quantization (in this sense dis- 
order is strong for them), we expect that their dynam- 
ical response has only weak interaction corrections (for 
kpas ^ 1 and kpRc ^ 1). On the experimental side, 
the CR of strongly correlated particles has been inves- 
tigated in the high-B limit (e.g., ||3l|-|3^ and references 
therein], but not for the translational-symmetry broken 
state in weaker magnetic fields. 

Secondly, because of electron-electron interactions, the 



quantity that is directly measured in far- infrared (e.g., 
psf ) or microwave (e.g., |3^-|3^]) experiments may be re- 
lated to the conductivity in a complex way. In the pa- 
per, we have provided explicit results for the conductiv- 
ity; more specifically, for Sciu>) oc Keaxxi^j) +hnaxy{LLi)- 
We expect only small interaction corrections to Sc{i^)- 
It is important, however, that the conductivity expresses 
the current as a response to the total (screened) elec- 
tric field. As such, it is given by the irreducible, with 
respect to the Coulomb interaction, density-density re- 
sponse function K{q,q',uj). On the other hand, in con- 
trast to dc measurements, what is probed directly in the 
ac transmission experiments is a response to the exter- 
nal field. The latter is given, in an operator form, by 
the reducible polarization e~^K, where e is the dielectric 
function. Put differently, the dissipated power is mea- 
sured in units of the intensity of the incident wave. Since 
the measured absorption is increased near zeros of e, i.e., 
on resonance with plasma oscillations, the dynamical re- 
sponse of a 2DEG calculated in the paper is in general 
masked in the absorption experiments by the excitation 
of magnetoplasmons (for optical experiments with mag- 
netoplasmons see [ ^5p9| - ^ , ^ , ^3| -^ ) . The measured 
quantity appears to depend in an essential way on the 
experimental setup (which may be very different, cf. ]3^ ] 
and ||,||). 

Let us mention one more point related to edge magne- 
toplasmons in AD arrays, i.e., a soft mode of plasma os- 
cillations localized near the sharp edges of ADs ( p6| , ^ , 
for experiments see [|9|-|2]j2^,^-|5[ ) . It is often asserted 
that there is an intimate connection between the collec- 
tive edge excitations and skipping orbits. In fact, this 
notion may be very misleading. Indeed, the width of the 
strip around an AD in which the skipping orbits propa- 
gate along the edge is 2i?c. However, the current in the 
edge magnetoplasmon mode decays, in general, on a dif- 
ferent scale: e.g., at ^ Wc and Rc 3> as this scale is 
given by R^/ub ^ Rc- Hence, the main contribu- 
tion to the edge magnetoplasmon current may come from 
electrons that do not at all collide with the AD. Since the 
frequency dependence of the dissipative response of the 
skipping orbits and that of the edge magnetoplasmons 
belong to different ranges of w, it appears to be quite 
possible that they can be measured separately in trans- 
mission experiments. 

In conclusion, we have seen that the dynamical con- 
ductivity of electrons scattered on impenetrable ADs in 
the presence of a magnetic field reveals strong memory ef- 
fects in the electron dynamics, associated with skipping 
orbits bound to ADs. These lead to the CR lineshape 
which is not at all characterized by the Drude scattering 
rate. The contribution of the skipping orbits Sc{uj) is 
broadened on a scale of the cyclotron frequency lOc and 
vanishes at lOc in a nonanalytical way as |w — Wc|- Apart 
from these two features, Sc{uj) exhibits different behav- 
ior depending on the ratio of the cyclotron radius Rc and 
the AD radius a. At large Rc/a, Sc{(^) oscillates with a 
period lUc up to w = uJcRc/ a and shows a series of square- 
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root spikes for larger w. At small Rc/a, Sdw) has a hard 
gap between two sharp peaks located at w ~ oJcRc/cl and 
Lu = ujc- We hope that these results will stimulate further 
experimental work on the ac conductivity in AD arrays. 



ACKNOWLEDGMENTS 

We thank LW. Engel, R. Haug, F. Hohls, I.V. 
Kukushkin, U. Merkt, R.J. Nicholas, and J.H. Smet for 
discussions concerning the experiments. This work was 
supported by SFB 195 and the Schwerpunktprogramm 
"Quanten-Hall-Systeme" of the Deutsche Forschungsge- 
meinschaft, by INTAS Grants No. 99-1070 and 99-1705, 
and by RFBR Grants No. 99-02-17093 and 00-02-17002. 



[1] 
[2] 

[3] 
[4] 
[5] 

[6] 
[7] 

[8] 
[9] 
[10] 

[11] 
[12] 
[13] 
[14] 

[15] 



Also at A.F. loffe Physico- Technical Institute, 194021 
St. Petersburg, Russia. 

D.K. Ferry and S.M. Goodnick, Transport in Nanostruc- 
tures (Cambridge University, Cambridge, 1997). 

C. W.J. Beenakker and H. van Houten, in Semiconductor 
Hetero structures and Nano structures, Solid State Phys. 
44, 1 (1991). 

D. Weiss, G. Liitjering, and K. Richter, Chaos, Sohtons 
and Fractals, 8, 1337 (1997). 

R. Fleischmann, T. Geisel, R. Ketzmerick, and 
G. Petschel, Physica D 86, 171 (1995). 
G.M. Gusev, P. Basmaji, Z.D. Kvon, L.V. Litvin, 
Yu.V. Nastaushev, and A.I. Toropov, J. Phys.: Condens. 
Matter 6, 73 (1994). 

K. Tsukagoshi, S. Wakayama, K. Oto, S. Takaoka, 
K. Murase, and K. Game, Phys. Rev. B 52, 8344 (1995). 
G. Nachtwei, Z.H. Liu, G. Liitjering, R.R. Gerhardts, 
D. Weiss, K. von Klitzing, and K. Eberl, Phys. Rev. B 
57, 9937 (1998). 

O. Yevtushenko, G. Liitjering, D. Weiss, and K. Richter, 
Phys. Rev. Lett. 84, 542 (2000). 

A.D. Mirlin, J. Wilke, F. Evers, D.G. Polyakov, and 

P. Wolfle, Phys. Rev. Lett. 83, 2801 (1999). 

A.D. Mirlin, D.G. Polyakov, F. Evers, and P. Wolfle, 

cond-mat/0104096[ 

D.G. Polyakov, F. Evers, A.D. Mirlin, and P. Wolfle, 



cond-mat/0104417 



M.M. Fogler, A.Yu. Dobin, V.I. Perel, and B.I. 
Shklovsku, Phys. Rev. B 56, 6823 (1997). 
F. Evers, A.D. Mirlin, D.G. Polyakov, and P. Wolfle, 
Phys. Rev. B 60, 8951 (1999). 

E.M. Baskin, L.N. MagariU, and M.V. Entin, Zh. Eksp. 
Teor. Fiz. 75, 723 (1978) [Sov. Phys. JETP 48, 365 
(1978)]; E.M. Baskin and M.V. Entin, Physica B 249- 
251, 805 (1998). 

A.V. Bobylev, F.A. Maa0, A. Hansen, and E.H. Hauge, 
Phys. Rev. Lett. 75, 197 (1995); J. Stat. Phys. 87, 1205 
(1997). 



[16 
[17 

[is: 

[19 

[20 
[21 

[22; 

[23; 

[24 

[25' 
[26; 

[27 

[28; 

[29 

[30 

[31 

[32; 

[33 

[34; 



[35] 
[36] 
[37] 
[38] 



A. Kuzmany and H. Spohn, Phys. Rev. E 57, 5544 

(1998) . 

M.M. Fogler and B.I. Shklovskii, Phys. Rev. Lett. 80, 
4749 (1998). 

T. Ando, J. Phys. Soc. Jpn. 38, 989 (1975). 

D. Weiss, M.L. Roukes, A. Menschig, P. Grambow, 
K. von Klitzing, and G. Weimann, Phys. Rev. Lett. 66, 
2790 (1991). 

R. Fleischmann, T. Geisel, and R. Ketzmerick, Phys. 
Rev. Lett. 68, 1367 (1992). 

E. Vasiliadou, R. Fleischmann, D. Weiss, D. Heitmann, 
K. von Klitzing, T. Geisel, R. Bergmann, H. Schweizer, 
and C.T. Foxon, Phys. Rev. B 52, R8658 (1995). 

J. Wilke, A.D. Mirlin, D.G. Polyakov, F. Evers, and 
P. Wolfle, Phys. Rev. B 61, 13774 (2000). 
M.M. Fogler, A.A. Kulakov, and B.I. Shklovskii, Phys. 
Rev. B 54, 1853 (1996). 

R. Moessner and J.T. Chalker, Phys. Rev. B 54, 5006 
(1996). 

W. Kohn, Phys. Rev. 123, 1242 (1961). 

H. Fukuyama, Y. Kuramoto, and P.M. Platzman, Phys. 

Rev. B 19, 4980 (1979). 

C. Kallin and B.I. Halperin, Phys. Rev. B 31, 3635 
(1985). 

E. Fradkin and S.A. Kivelson, Phys. Rev. B 59, 8065 

(1999) . 

A. H. MacDonald and M.P.A. Fisher, Phys. Rev. B 61, 
5724 (2000). 

F. von Oppen, B.I. Halperin, and A. Stern, Phys. Rev. 
Lett. 84, 2937 (2000). 

B. A. Wilson, S.J. Allen, Jr., and D.G. Tsui, Phys. Rev. 
B 24, 5887 (1981). 

M.J. Chou, D.C. Tsui, and G. Weimann, Phys. Rev. B 
37, 848 (1988). 

M. Besson, F. Gornik, CM. En- 

gelhardt, and G. Weimann, Semicond. Sci. Technol. 7, 
1274 (1992). 

G. M. Summers, 
R.J. Warburton, J.G. Michels, R.J. Nicholas, J.J. Har- 
ris, and C.T. Foxon, Phys. Rev. Lett. 70, 2150 (1993). 
Note that a weak dependence of the cyclotron frequency 
on spin orientation may become important for the ex- 
planation of the experimentally obtained lineshape, see 
N.R. Cooper and J.T. Chalker, Phys. Rev. Lett. 72, 2057 
(1994). 

E. Batke, D. Heitmann, and C.W. Tu, Phys. Rev. B 34, 
1986. 

L.W. Engel, D. Shahar, C. Kurdak, and D.C. Tsui, Phys. 
Rev. Lett. 71, 2638 (1993). 

F. Hohls, U. Zeitler, and R. Haug, Phys. Rev. Lett. 86, 
5124 (2001). 

P.D. Ye, L.W. Engel, D.C. Tsui, J.A. Sim mons , 
J.R. Wendt, G.A. Vawter, and J.L. Reno, cond 
mat/0103127l 



[39] K. Kern, D. Heitmann, P. Grambow, Y.H. Zhang, and 
K. Ploog, Phys. Rev. lett. 66, 1618 (1991). 

[40] Y. Zhao, D.C. Tsui, M. Santos, M. Shayegan, R.A. Ghan- 
bari, D.A. Antoniadis, and H.I. Smith, Appl. Phys. Lett. 
60, 1510 (1992). 

[41] A. Lorke, I. Jejina, J. P. Kotthaus, Phys. Rev. B 46, 



8 



12845 (1992). 

[42] K. BoUwcg, T. Kurth, D. Hcitmanri, E. Vasiliadou, 

K. EbcrL and H. Brugger, Phys. Rev. B 52, 8379 (1995). 
[43] K. Bollweg, T. Kurth, D. Heitmann, V. Gudmundsson, 

E. Vasiliadou, P. Grambow, and K. Eberl, Phys. Rev. 

Lett. 76, 2774 (1996). 
[44] M. Hochgrafe, R. Krahne, Ch. Heyn, and D. Heitmann, 

Phys. Rev. B 60, 10680 (1999). 



[45] S. Cina, D.D. Arnone, H.P. Hughes, C.L. Foden, 
D.M. Whittaker, M. Pepper, and D.A. Ritchie, Phys. 
Rev. B 60, 7780 (1999). 

[46] V. Fessatidis, H.L. Cui, and O. Kiihn, Phys. Rev. B 47, 
6598 (1993). 

[47] S.A. Mikhailov and V.A. Volkov, Phys. Rev. B 52, 17260 
(1995); S.A. Mikhailov, ibid. 54, 14293 (1996). 



9 



